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Abstract 

It is argued that the discovery of cosmic acceleration could have been anticipated on thermody- 
namic grounds, namely, the generalized second law and the approach to equilibrium at large scale 
factor. Therefore, the existence of dark energy -or equivalently, some modified gravity theory- 
should have been expected. In general, cosmological models that satisfy the above criteria show 
compatibility with observational data. 
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I. INTRODUCTION 

The standard cold dark matter (SCDM) model [l| was in good health until around the 
last decade of the previous century when it became apparent that the fractional density of 
matter falls well below the Einstein-de Sitter value, Vt m = 1 -see e.g. 

as 

. The death blow 

came at the close of the century with the discovery of the current cosmic acceleration 
something the said model cannot accommodate by any means. However, to account for the 
acceleration in homogeneous and isotropic models one must either introduce some exotic 
energy component with a huge negative pressure (dubbed dark energy) or, more drastically, 
devise some theory of gravity more general than Einstein relativity |5j. Thus, both solutions 
appear somewhat forced and not very aesthetical. Here we argue that dark energy (or 
something equivalent) is demanded on thermodynamic grounds. In other words, we provide 
what we believe is a sound thermodynamic motivation for the existence of dark energy. 

Our argument is based on that the natural tendency of systems to evolve toward ther- 
modynamical equilibrium is characterized by two properties of its entropy function, S(x), 
namely, it never decreases, dS(x)/dx > 0, and is convex, d 2 S(x)/dx 2 < |6j. In the context 
of an ever expanding Friedmann-Roberson- Walker (FRW) cosmology this translates in that 
the entropy of the apparent horizon plus that of matter and fields enclosed by it must fulfill 
S'(a) > at any a -the generalized second law (GSL)- as well as S"(a) < as a — > oo, where 
a is the scale factor of the FRW metric and the prime means d/da. The apparent horizon 
in FRW universes always exists (which is not generally true for the particle horizon and the 
'uture event horizon) and is known to posses not only an entropy proportional to its area 



boundary [10 ]. 



7|, |8[ but also a temperature js]. Besides, it appears to be the appropriate thermodynamic 



Before we proceed, it is fair to recall the (at least theoretical) existence of systems lacking 
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12] . The latter system 



any global maximum entropy state, such as Antonov's sphere 
consists in a sphere enclosing a number of particles that share some total energy. If the sphere 
radius happens to increase beyond some critical value, the system becomes unstable and the 
entropy function ceases to have a global maximum. We shall apply our argument under the 
assumption that the Universe tends to a state of maximum entropy irrespective of whether 
it may eventually reach it or not. It would be odd and frustrating that such a universal 
principle as the second law of thermodynamics could not be applied to the Universe as a 



whole, especially given the close connection between thermodynamics and gravity [13J, Il4j . 

Section II illustrates why dark energy (or some or other modified gravity model) is re- 
quired on thermodynamic basis and study some dark energy models to see whether they 
fulfill the thermodynamic criteria. Section III applies the said criteria to some representative 
modified gravity models. Finally, section IV summarizes and discusses our findings. As is 
customary, a naught subscript stands for the present value of the corresponding quantity. 



II. WHY DARK ENERGY WAS TO BE EXPECTED 

Let us consider a FRW universe. Its entropy is contributed by two terms: the entropy 
of the apparent horizon which is proportional to its area, A = 4:irr\, and the entropy of 
the fluids enclosed by the horizon. Here ta = (\/H 2 + {k/a 2 ))~ l denotes the radius of the 
horizon and H the Hubble factor of the FRW metric [7]. As is well-known, Sa = ^f- -gr 
where £pi and ks stand for Planck's length and Boltzmann's constant, respectively. 

In virtue of the first Friedmann equation 



k 



3H 2 +3— = 8vrGp, (1) 
a 2 



where p is the total energy density, we can write 



as well as 



= ( 2) 



A'- — 1 +w (3) 
A ~2G ap ' {6) 



where we have used the conservation equation p' = —3(1 + w)p/a, where w = p/p denotes 
the overall equation of state parameter, i.e., not just of dark energy. From §5§ the area will 
augment in expanding universes if 1 + w > and decrease otherwise. 
A further derivation with w = constant yields 

x = -|^W (ap ' + pH 5^v (1 + 1 " ,(2 + 3t " , • (4) 

Accordingly, A" < for — 1 < w < —2/3, and A" > otherwise. Thus, the above criterion 
disfavor the dominance of fluids at late times such that the overall equation of state is either 
of phantom type or larger than —2/3. 



We note in passing that when w ^ constant last equation generalizes to 



A 
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w + 



(1 + w) (2 + 3w) 



(5) 



2Gap 



a 



Let us consider the entropy associated to the fluid enclosed by the apparent horizon. If 
the fluid is just dust (pressureless matter, subscript m), i.e. p = p m and p = p m = 0, the 
unphysical result that the temperature vanishes follows whence Euler's equation, T s = p + p, 
cannot be used to determine the entropy of cold matter (i.e., dust). We proceed instead 
as follows. Every dust particle contributes a given bit -say ks- to the fluid entropy. So, 
within the apparent horizon we will have S m = ks N being N = (4ii/3)r\n the number of 
particles there, and n = Uq a -3 the (conserved) number density of dust particles. Then, 



enclosed by the horizon because the latter encompasses more and more particles as the 
universe expands (i.e., as H diminishes). Thus, S'^ + S'^ > 0. 

We may now readily understand why sooner or later the Universe should accelerate, i.e., 
why it must be endowed with dark energy (subscript x, equation of state — 1 < w x < —2/3), 
or something dynamically equivalent at the background level (as a suitably modified gravity). 
We note parenthetically that, at late times, dark energy is to dominate over all other energy 
components whereby w ~ w x as a — > oo. Were the Universe dominated by radiation and/or 
matter for ever, S'^ could never become negative. And if the horizon entropy dominated the 
total entropy, then the Universe would never tend to a state of maximum entropy (compatible 
with the constraints of the system). 

Figure Q] sketches the evolution of the horizon area for a spatially flat FRW universe 
dominated by a mixture of cold matter and dark energy, i.e., A oc [f2 m0 a~ 3 + Vl x0 a~ 3( > 1+Wx * ) ]~ 1 
for w x = —5/6 (solid line), w x = —1 (dashed), and w x = —1.2 (dot-dashed) -in the three 
cases we have assumed Q m0 = 0.3. It is readily seen that for —1 < w x < —2/3 (solid 
and dashed lines) the curvature evolves from positive to negative values and stays negative 
forever. By contrast, when the dark energy component is phantom, w x < —1 (dot-dashed 
line) the curvature goes from positive to positive values with an intermediate, transient, 
phase of negative values about the maximum. Clearly, while in the two first instances the 
Universe tends to thermodynamic equilibrium at late times in the case of phantom it does 




Hence oc > 0. We note parenthetically that the entropy augments in the volume 
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not. Needless to say, in the absence of dark energy (not shown) we would have A" > for 
whatever value of the scale factor. 
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Figure 1. Schematic evolution of the area of the apparent horizon in a FRW universe dominated 
by cold matter and dark energy with w x = —5/6 (solid line), w x = — 1 (dashed), and w x = —1.2 
(dot-dashed). In plotting the graphs we have taken Q m Q = 0.3. 

For completeness, Fig. [2] illustrates the influence of the spatial curvature on the evolution 
of A when a ^> 1 for three ACDM models (w x = — 1). As can be seen, the curvature of the 
graph for positively curved models remains positive (recall that Qk = —k/(a 2 H 2 )). 

We believe that the GSL alongside the criterion that the total entropy S must fulfill 
5"' < at late times may serve to discard some cosmological models and to set limits on the 
evolution of w x . 

It will prove useful to express A and its derivatives in terms of the deceleration parameter, 
q = —d/(aH 2 ). The latter can be written as 




(7) 
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Figure 2. Schematic evolution of the area of the apparent horizon in a FRW universe dominated 
by cold matter and cosmological constant (uu x = — 1). From top to bottom, f2&o = —0.02, 0, and 
0.009, respectively, -see [jj]]. For a values of the order unity and below (not shown) the graphs 
practically overlap each other. In plotting the graphs we have chosen Qao = 0.7. 



In spatially flat FRW universes A = 4irH 2 , whence A' = 8n(l + q)/H 2 and 



A" = 2A 



;i + q)(l + 2q) + 



(8) 



It follows that for q > —1/2 the second derivative of the horizon area, A", cannot go from 
positive to negative values unless q' < 0. This excludes evolutions of q as in Fig. |3j The 
latter corresponds to cosmological models in which the current accelerated expansion is just 
transitory -after a period of dark energy dominance, pressureless matter is assumed to take 



over 
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At this stage it seems fitting to ask ourselves whether the second derivative of the entropy 
of the energy components -that for simplicity we will call fluids, subscript /- enclosed within 
the apparent horizon will be positive enough so that the sum S'j + S'^ be positive. However, 
before proceeding we wish to remark that if the energy component is a scalar field in a pure 
quantum state, or the cosmological constant, it will have no entropy at all. Nevertheless, we 
wish to consider the possibility of the scalar field being in a mixture state and accordingly 
entitled to have entropy (obviously, this excludes the case Wf = —1). 

For simplicity we will take k = and a single fluid component with Wf = constant. The 
entropy of the fluid filling the volume enclosed by the apparent horizon can be estimated by 




a 



Figure 3. Schematic of the evolution of the deceleration parameter, q(a), given by Eq. ©, in 
a FRW universe where the present acceleration stage is transitory and reverts for ever to matter 
domination. 



virtue of Gibbs equation, 



47T 



T f dS f = d[p fY H- 

With help of and noting that p = p/ 

dS f _ 



2GT 



f 



da 



'47T 

u>fPf d [ y H ~ 3 



1 M o 



(9) 



Using the second Friedmann equation. 

dH 



da 



-AnG (1 + w f ) 



Pf 
aH 



(10) 



and Eq. ([T]) once more, we get 

7GT f 



dSf 3 . . . . 1 

--(l + Wf )(l + 3w f )—. 



da 2 

Thus, dSf/da will be negative or nil for — 1 < Wf < — 1/3 and positive otherwise. 
From last equation we can write, 

d 2 S f 3 , , , v d 



(11) 



2G 



V 
da 2 



-(l +Wf )(l + Z Wf )-(aHT f )- 1 . 



The evolution of the temperature of a perfect fluid, dlnTf /dlna = —3vjf, readily follows 
from Gibbs' equation and the condition for dSf to be a differential expression -see e.g. 
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181 ] . Recalling that Wf = constant, it leads to Tf = Tf Q a 3w f. Using now pf 



Pfo 



a 3 ( 1+M '/) we arrive to 



2G 1^ = iT^Ho (1 + Wf) (1 + 3Wf) (1 + 9Wf) a{9Wf ~ m ■ (12) 

For — 1 < Wf < —1/3 (dark energy) as well as for Wf > —1/9 one follows that S'f > 0. 
However, S'f decreases with expansion for Wf < 1/9. The question arises whether the sum 
S'f + S' A will be positive or negative in the long run. To answer it we inspect the ratio 
between the right hand sides of Eqs. ( !T2|) and (@| and note that w = Wf, 

S oc G 3K-D/2 . (13) 

This expression vanishes for Wf < 1 as a — > oo, whereby S'f + S'^ < in the long run 
provided — 1 < Wf < — 2/3, i.e., if the fluid is dark energy. 

It remains to be seen whether the GSL, S'f + S' A > 0, is fulfilled. Note that in virtue of 
Eqs. © and (dU), with T f oc a~ 3w f, the ratio 

|f oc 1/(flT/g) oc pf oc a 3 ^- 1 )/ 2 (14) 
S' A 1/apf 1 

vanishes in the long run provided Wf < 1. Thus albeit S'f is negative for —1 < Wf < —1/3, 
the GSL is satisfied for Wf > —1 when a — > oo. 

Altogether, dark energy with constant equation of state in the interval (—1, —2/3) respects 
the GSL as well as the criterion that S'f + S' A < when a — > oo. This result is consistent 
with the tightest observational constraints available -cfr. Table in IV in In our view, 
cosmological models that meet both criteria should be preferred to those failing any of the 
two. 



A. Parameterized w x models 

In this subsection we consider the model of Barboza and Alcaniz in which the dark energy 
equation of state is parameterized in terms of redshift as 2(]] 

/ \ -2(1 + Z) 

w x {z) = W + W!—— — — . (15) 
1 + z z 

It has the advantage over other parameterizations of not diverging at any time. We recall, 
parenthetically, that the redshift is related to the scale factor byl + z — 1/a. 
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The pair of free, constants parameters wo and W\ in ( fT5|) ought to be restricted by physical 
requirements and observation. A first constraint wq + w\ < follows from demanding that 
dark energy be subdominant at early times (when z> 1), otherwise cosmic structure would 
never had formed. Next, we will use the GSL together with the condition that the sum 
S'm + S" + S'a. be negative or nil to set further constraints on the said pair. 

We begin by writing the Hubble function of a spatially-flat FRW universe dominated by 
cold matter and dark energy, with equation of state f|T5|) . 



H 2 = Hi [Q m0 a- 3 + Q x0 a-^ 1+W0+Wl \2a 2 - 2a + if^l 2 ] . (16) 

We next consider the entropy of the apparent horizon, Sa oc A. Two separate cases arise: 
(i) wi < 2/3 and —2/3 < w < —w\ as well as (ii) wo < —2/3 and wo + w\ < 0. In the first 
case, S' A grows without bound and we do not consider it any further. In the second one S' A 
approaches zero asymptotically from above. As a — > oo, we see that w x (a) — > wo, fl x — > 1, 
aw' x — > (~ a -1 ) and H 2 ~ a _3< - Wo+1 ^ whence all the terms but the last one, that decays as 
a -1 , in the second derivative of the horizon entropy, 



rill All 1 

b A OC A CX 



53 i 

3(w x £l x ) 2 + ~w x fi x - -w 2 x Q x + 1 + -a Q x w' T 



(17) 



a 2 H 2 

present a like behavior. In consequence, the quantity in the square brackets tends to a 
constant, and S A ~ a 3wo+l as a — > oo. For S A to be negative that constant must be negative 
as well, then 3wq + 5w + 2 < 0. 

Finally, this constraint and the previous one set the accessible range to 

Wx < 2/3 f 2/3 < wi < 1 

' U I ' (18) 

-1 < w < -2/3 1 -1 < w < -wx, 
as shown in Fig. (j4j). 

Obviously, when the entropy of matter and dark energy are taken into account the pre- 
vious result may vary. With the help of the equation for the temperature of dark energy 

T x {a) = T x0 a- 3{wo+Wl \2a 2 - 2a + l) 3wi/2 , (19) 

it can be checked that 

S' S" 

~ d and -f- ~ C 2 , 

where C\ and C*2 are constants. Hence we focus on matter contribution with respect to the 
horizon. 
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Figure 4. The available region in the (wq,wi) plane for the parameterized-u> x model of Barboza 



and Alcaniz 



20|. 



When w > it follows that 



qi q// 

°m °m „-3/2 

Q/ Q" ' 



whereby the horizon contribution prevails in the long run. However, we have already dis- 
carded this solution because S' A diverges. In a sense, this is a consistency check on our 
criteria because observationally w7o is known to be negative. 
In the opposite case, wo < 0, one obtains 

q/ q// 

°m °m 3(mj -1)/2 
ci on 

Again, the horizon contribution prevails, and the result f lT8|) stands. 

The FRW universe model dominated by cold matter plus dark energy, with equation of 
state parameterized according to Eq. (1151) . has been observationally constrained in 20[ by 
using data from supernovae type la, baryon acoustic oscillations, and cosmic background 
radiation. The best fit values of the parameters at la confidence level were found to lie in 
the ranges 

- 1.35 < w < -0.86, -0.33 < wi < 0.91, (20) 



in full consistency with the above result, Eq. ([18]) . 
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B. Chaplygin gas model 



The original Chaplygin gas model unifies matter and dark energy in the sense that they 
are no longer two separate components but a unique entity that mimics cold matter at early 
times and a cosmological constant at late times 21] . Its equation of state p = —A/p, with 
A a positive constant and p the energy density of the gas (i.e., the total energy density), is 
obtainable from the Nambu-Goto action for a d-brane moving in the d + 1 dimensional bulk 



221 ] . In a FRW universe, the dependence of the energy density on the scale factor reads 



p = \J A + (B/a 6 ), where B is nonnegative integration constant; or equivalently, 

p = Po [l - £1* + ft* a~ 6 ] 1/2 (21) 

with = B/(A + B). 

The evolution of the horizon area (A oc 1/p) is akin to the one depicted by the solid 
line in Fig. [T| that is to say, it grows monotonously and presents negative curvature at 
sufficiently large scale factor. 

By contrast, since the entropy of the Chaplygin gas obeys Sen oc (1 + w)/(HT), where 

(22) 



p 

w — — 
P 



■ a. 



1 - Si, 

refers to the equation of state parameter of the gas (which evolves from zero to —1), and 
the temperature is governed by T'/T = — 3w/a, one follows that S' Ch < and S 1 ^ > for 
large scale factor. However, as it can be checked, S' Ch / ' S' A — > and S'^ h /S' A — > as a — > oo. 
Accordingly, the total entropy, S = Sen + Sa, fulfills S' > and S" < in the same limit. 
In other words, the GSL is respected in that limit and the Universe tends to thermodynamic 
equilibrium in the long run. 



C. Holographic, interacting models 

Here we first consider a spatially-flat, holographic, interacting model dominated by pres- 
sureless matter and dark energy. The latter component is assumed holographic in the sense 
that its energy density varies as the area of the apparent horizon and interacts with matter 
at a given, non-constant, rate. As a result w x decreases with expansion and the fractional 
densities of both components remain constant. This much alleviates the cosmic coincidence 



problem |23[ and the model shows compatibility with observation 
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24] 



Inspection of the Hubble factor, given by Eq. (2.4) of 24 ] 

r /. r 



Ha) = H n 



+ 1 - 



-3/2 



(23) 



3H r V. 3# r, 

readily reveals that the graph of the evolution of the area of the apparent horizon, A oc 1/H 2 , 
has positive curvature at the beginning of the expansion and gently evolves to negative 
curvature values remaining thus for ever, similarly to the solid line in Fig. [1], for the best fit 
observational values, T / H = 0.563 and r = 0.452, obtained in [24]. 

As can be readily checked, the second derivative of the entropy of the dark energy within 
the horizon, S x , is negative as well. Effectively, 

l + w x 



4tt„ 3 _ 4tt 1 l + w x 

x 3 rAS 3 m n 



Px OC 



(24) 



HT X ' 

where r = p m /p x and T x oc a~ 3v)x when a — > oo. For a ^> 1 one has that w = constant, and 
one obtains S x cx 3(1 + w x ) w x (3w x — 1) < since, in the model at hand, 1 + w x < 0. 
As for the entropy of pressureless matter 



On 



l 4?r ~3 

k- 



r a n 



47T p r , 



k 



B 



AlTx 1 r 1 
3 mii^ 1 + r i/ 



3 " 3 mi/ 3 

where m stands for the mass of the particles. (Recall that in this model r = constant). 
We then have, S' m oc -H'/H 2 and S'^ oc 2(H' 2 /H 3 ) - (H"/H 2 ). Using (EHJ we find that 
S 1 ^ > and 

S" oc 7 /3 2 



, 2 7 - 5 a 3 / 2 
a 5 ii 3 



where /3 is a short-hand for the dimensionless combination r/(3iio r ) and 7 = (1 — /3)//3; 
both quantities lie in the range (0,1). For a ^> 1, the second term is the leading one 
whence, at sufficiently large scale factor, S'^ is also negative. In summary, in this model, 
S'a + S'i + S'x 1 < for a -> oo. 

It would be misleading, however, to believe that for all holographic models A" < when 



a ^> 1. For instance, in the holographic model of Gao et al. 
length as infrared cutoff, the Hubble equation reads 



25] , which uses the Ricci's 



H = H x Q k o a 2 + Q m0 a 3 + fi r0 a 4 + 



(26) 



2 - ft 

where the subscript r stands for radiation; a ~ 0.46 and /o ~ 0.65 are dimensionless 
parameters. As Fig. |5] shows, the area of the apparent horizon grows to a maximum to 
monotonously decrease afterwards for ever. Consequently, A" > as a — > oo. 
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Figure 5. Schematic evolution of the area of the apparent horizon of the holographic, non- 



interacting model of Gao et al. 



251 ] . The curvature is negative about the maximum and becomes 



larger than zero later on and remains thus for ever. In drawing the graph we have adopted the 
values in the said reference: Qko = 0, U r o = 8.1 x 1CT 5 , Q m o = 0.27, Q x o = 0.73, a = 0.46, and 
/o = 0.65. 

On the other hand, the entropies of the fluid components (radiation, matter, and dark 
energy) decrease with expansion and their second derivatives are positive for large scale fac- 
tor. Thus, this model violates the GSL and does not approach thermodynamic equilibrium 
at late times. 

III. MODIFIED GRAVITY MODELS 



Models that depart from Einstein gravity may lead to a late acceleration era without 
the help of any exotic component of negative pressure. Here we examine some of them, 
namely: the one based in the brane-induced gravity model of Dvali et al. 26] (see 27], |28|). 
the original Cardassian model proposed by Freese and Lewis [29f| , and the torsion model of 
Bengochea and Ferraro 30|, to check whether they fulfill these criteria. Before we proceed, 
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it is fair to say that, as far as we know, it has not been rigorously proven that the entropy 
associated to the apparent horizon of the two latter models simply obeys Sa oc A, but it 
seems to us a very reasonable assumption and we shall adopt it. 



Dvali-Gabadazze-Porrati's model 



This model considers our 4-dimensional Universe as a brane embedded in a 5-dimensional 
bulk with flat Minkowski metric. As a consequence of the brane-induced term, the conven- 
tional Friedmann's equation modifies to 

1 



H z + — 



P 



3M* pl 



(27) 



where p and r c stand for the total energy density (matter plus radiation in this model) and 
the crossover scale below which gravity appears as four dimensional, respectively. 



The entropy of the apparent horizon, as computed by Sheykhi et al. 



31] 



s A = k B 3wf * 



02 



1 + ^ 

TV, 



(28) 



is no longer proportional to the horizon area but contributed by two terms. The first one 
corresponds to the Bekenstein-Hawking entropy and accounts for gravity on the brane. The 
second one is related to the bulk. Its evolution in terms of the scale factor is sketched in 
the left panel of Fig. |6j From a certain point on it increases monotonously with negative 
curvature. As before, the entropy of dust and radiation are proportional to (aH)~ 3 . Both 
of them decrease monotonously with positive curvature. It is found that 

S' S' 1 



— )• 



Ho 



(29) 



(subscript r for radiation), bear in mind that r c has physical dimensions of time. Thus, in 
the long run it leads to a negative constant -though the latter does not exactly coincides in 
the radiation and matter cases. 

Once again the question arises whether the total entropy, S = S m + S r + Sa, obeys 
the GSL and presents negative curvature when a>l. To answer this consider the limiting 
expressions when a — » oo, 



Qf Cff 

*_> O A 



S" 



&A 



1 _l_ m I r 
° A °A. 

1 nil r~in 



Qll 



Qll 



-> 



-)> 0. 



(30) 
(31) 



14 



Figure 6. Left panel: Schematic evolution of the entropy of the horizon of the Dvali-Gabadazze- 
Porrati model. Right panel: The single graph depicts the qualitative evolution of the entropy of 
radiation as well as of cold matter inside the horizon. In drawing the graphs we have used the 
best-fit values of the parameters of the model, t c Hq = 1.21 and O m o = 0.18, obtained in 32]. 



The first one tends to zero from above, the second one from below (following the healthy 
behavior of the apparent horizon entropy). Both relations when taken together imply the 
inequality 



1 



k B n + 4 



PrO c 



rO 



(32) 



which translates into an upper bound on the current number density of dust particles, 



n < 



4^ 



10 38 cm -3 



(33) 



Since it is fulfilled by a huge margin the GSL is satisfied and S" results negative in the long 
run. 



B. Cardassian model 

In this spatially-flat FRW model the first Friedmann equation acquires an extra term 
that accounts for acceleration at sufficiently high redshifts, 

H*= 8 -^p + Bp*. (34) 

Here two new non- negative constants, B and a, appear while p stands for the energy density 
of cold matter -the only energy component. For a < 2/3 the universe features a transition 
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from q > to q < at redshift 



~1r 



1/3(1-") 



L(4ttG/3) P ^. 

without need whatsoever of dark energy. 

As can be checked, A = 1/ H 2 evolves in such a way that while A' > at all redshifts, 
A" changes from positive to negative values at some point and stays thus for ever. This was 
to be expected since the extra term on the right had side of ( |34|) dynamically amounts to 
the presence of some dark energy component. Thus, the overall behavior of A should be 
qualitatively similar to that of a model dominated by a mixture of pressureless matter and 
dark energy. To be more specific, comparison of the right hand sides of flM|) and the first 
Friedmann equation for a mixture of cold dark matter and dark energy with constant w x 
and k = 0, namely, H 2 = (8irG/3)[p m oa~ 3 + p x o a~ 3 ( 1+Wx '], shows that, at the background 
level, every Cardassian model can be mapped to a spatially-flat dark energy one satisfying 
a = 1 + w x and B = (8iiG/3) (Pxo/Pmo)- As a consequence A" < for sensible a values 
(i.e., < a < 2/3). 

Although the entropy of the matter enclosed by the apparent horizon varies as (aH)~ 3 , 
whence S' m < and S'^ > for a > 1, the ratios S'JA' oc a 3 ^" 1 )/ 2 and S'JA" oc a 3 ^ 1 )/ 2 
tend to zero in the same limit. It implies that S' m + S' A > and S'^ + S A < as a — > oo. 



C. Torsion model 



The torsion model of Bengochea and Ferraro 30] falls into a class of models that describes 
gravitation in terms of the torsion scalar r rather than curvature, R. Thus, the action takes 
the form 

I = j d^^V^ ( r + f( r ))+ Imatter , (35) 

where /(r) is a free function to be constrained by observation and experiments. An impor- 
tant advantage of this set of theories is that the field equations are second order as opposed 
to the four order equations of f{R) gravity. However, up to now, specific models have been 
proposed in the cosmological context only. Works investigating the symmetries and dynam- 
ics of the theory show that it can exhibit extra degrees of freedom since the theory is not 
local Lorentz invariant |33i ] . 
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Bearing in mind that for the FRW metric r = —6H 2 , from (1351) the generalized Friedmann 
equations 

H 2 = 8 -^p - f - - 2H 2 f T , (36) 
dH 2 _ 16ttG P + 6H 2 + f + 12H 2 f T 

din a ~ 2AH 2 f TT -2-2f T ' ^ ' 

follow. Here, p and p denote the energy density and pressure of the fluid component -which, 
we will assume, just cold matter. 
In the model of Ref. 30j 

f(r) = -a(-r)- n , (38) 

where a = (1 - Vt m0 ){<oH 2 ) l - n / \2n - 1). 

We evaluate the entropy of the apparent horizon and matter entropy using the best fit 
values of the parameters, n = —0.1 and f2 m0 = 0.27 -cfr. Ref. [30J]. Left panel of Fig. [7]shows 
the evolution of area of the apparent horizon in terms of the scale factor. The plot starts 
with positive curvature and is ever increasing. Later on, the effect of the "dark torsion" 
is felt and the sign of the curvature changes to remain negative for ever. The right panel 
qualitatively depicts the evolution of the entropy of cold matter (proportional to (aH)~ 3 ). 
It begins increasing to decrease after the maximum. The curvature is negative about the 
maximum only. Thus, the question arises whether for a — > oo the GSL will be respected and 
the sum S'^ + S'^ will be negative. A numerical study shows that the ratios \S' m \/S' A and 
^to/I^aI increase with the scale factor when a ^> 1. Therefore the answer to both questions 
appears to be no. 

Before closing this subsection it is fair to recall that it is still unknown whether /(r) mod- 
els admit black holes solutions. If they don't, it will become unclear that some connection 
entropy-horizon area really holds in these theories. 



IV. DISCUSSION 



As we have argued, neither a radiation nor a cold matter dominated universe can tend to 
thermodynamic equilibrium in the long run. By contrast, dark energy dominated universes 
may; this holds true irrespective of whether the dark energy component has entropy or not. 
Accordingly, dark energy (or something dynamically equivalent at the background level, 
such as a suitably modified gravity theory) appears thermodynamically motivated. In other 
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Figure 7. Schematic evolutions of the area of the apparent horizon in the case of the f(r) model of 
Ref . [30] (left panel) , and the entropy of the cold matter enclosed by the apparent horizon (right 
panel). In plotting the graphs we employed the best fit values of the said model, n = —0.1 and 
f^mO = 0.27. 



words, any of these two kind of ingredients was to be expected on thermodynamic grounds. 
We, therefore, should not wonder that the Universe is accelerating. However, it does not 
mean that every accelerating universe is thermodynamically motivated; that is the case, 
for instance, of any phantom dominated expansion with w x = constant, and some modified 
gravity models. 

One may object that if our reasoning were valid, the Universe would have never ceased 
to innate as it would mean a transition from acceleration, A" < 0, to deceleration A" > 0. 
Clearly if the primordial inflation would have lasted for ever, big bang nucleo-synthesis would 
never have occurred, galaxies couldn't have come into existence, and so on -something in 
stark contrast with observation. However, this reasoning is rather incomplete as it leaves 
aside the huge entropy generated during the reheating process at the end of inflation. In 
this explosive and quasi- instantaneous event the inflaton field relinquishes all its energy in 
the form of matter and radiation and, as a consequence, the Universe sees its temperature 



enormously increased 34(. Since this huge amount of matter and radiation thermalizes (a 
necessary condition for primordial nucleosynthesis) both second derivatives, S'^ and S", are 
negative and, as a result, we may well have that S" + S'^ + S'^ < 0. Obviously, this will 
depend on the specific inflationary model and the particular reheating process involved, but 
we are not aware of any general argument against the fulfillment of this inequality. Specific 
calculations in this connection will be the subject of a future research. 
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Interestingly enough, cosmological models complying both with the GSL and the thermo- 
dynamic criterion that, at late times, the Universe should approach thermodynamic equi- 
librium appear to be compatible with observational data. Nevertheless, some models that 
do not comply with the said criteria (as some phantom models and some modified gravity 
models) seem also consistent with the said data. 

Our argument could be falsified if it were discovered that in the past (but after the 
radiation-dominated period set in) the Universe experienced one (or more) transitions from 
deceleration to acceleration and back. Cosmologies of the type have been proposed to explain 



the seemingly periodic distribution of galaxies with redshift -see, e.g. [35], [36j- and as an 
expedient to solve the cosmic coincidence problem [37J. However, recent studies on the 
impact of hypothetical transient periods of acceleration-deceleration on the matter growth 



and on the radiation power spectrum 



391 ] from the decoupling era, a ~ 10 5 , to a = 0.5 



practically discard these periods in the redshift intervals considered. 

If eventually it gets confirmed that the present phase of acceleration is to last forever, it 
may be seen as an indication that the Universe as a whole obeys the laws of thermodynamics 
(with the reservation that not all models that accelerate at late times comply with them). 
If, on the contrary, the Universe resumed a decelerated stage one should either call into 
question the validity of applying the said laws in a cosmological setting or wait for a later, 
and definitive, accelerating era. 

Altogether, it is for the reader to decide which possibility looks the less queer: dark energy 
(or modified gravity) or a universe that will never approach thermodynamic equilibrium. The 
remaining possibility, either dark energy or modified gravity combined with an increasing 
departure from equilibrium, involves two oddities rather than one. 
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